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ABSTRACT 


THE DYNAMIC PROPERTIES OF THE FULL JOURNAL BEARING, 
CONSIDERING CAVITATION 


by Richard A. Kaufer 


Submitted to the Department of Naval Architecture 
and Marine Engineering and the Department of Mech- 
anical Engineering on 20 May, 1966 in partial 
fulfillment of the requirements for the degree of 
Naval Engineer and the degree of Master of Science 
Inoucenanicalebngilncerinz, 


The unsteady motion of e circumferentially grooved journal 
bearing is analyzed and programmed on a digital computer, using 
the short bearing approximation of Reynolds! equation. The 
boundary conditions of continuity are used, and particular em- 
pnasis is placed on determining the effect of the size and 
shape of tne cavitation zone on the unsteady characteristics of 
the bearing. 


Tabuiated data consisting of bearing loads, load angles, 
ana spring and damping coefficients are presented for a wide 
range of bearing dimensions and operating conditions. A com- 
puter program is also written wnich can be used to obtain a 
wider range of values than are tabulated, and/or as a sub- 
routine in an analysis of a journal bearing supported system. 


The analysis shows that the dynamic properties of a jour- 
паі bearing are greatly influenced by the shape and extent of 
the cavitation zone, thus indicating a need for experimental 
work with the primary object of ascertaining film behavior 
under unsteady conditions. 


Thesis Supervisor: 


Herbert H. Richardson 


Associate Professor of Mechanical Engineering 
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NOMENCLATURE 


damping coefflcient, as in B,.; lb-sec/in. 
Bman ° 


dimensionless damping coefficient, 
pa LR 


radial bearing clearance; in. 
bearing center 

journal center 

perturbed journal center 
perturbation of journal center; in. 
journal diameter; in. 

journal eccentricity; in. 

total bearing load; lb. 


total steady state bearing load; lb. 
21 


pa LR 
bearing forces in 1 and 2 directions respectively; 15. 





dimensionless bearing, load, 


Steady state bearing forces in the 1 and 2 directions 
respectively; lb. 


bearing forces in r and s directions respectively; lb. 


radial bearing clearance at any point; c(l+ncos 8); 
in: 


radial bearing clearance along the upstream boundary 
of the cavítated zone; c(1+ncos0*); in. 


[nsa 


spring coefficient, as in K,.; lb/in. 


Бе 


p, LR 


excluding circumferential 


dimensionless spring coefficient, 


total bearing length, 
groove; in. 


eccentricity ratio, — 








х,у 


х,у,2 


pressure at any point; psia. 
feed pressure; psia. 
feed pressure, (Po ~ Pa): psig. 
ambient pressure; psia. 

Po 
pressure ratio, — 

а. 

volumetric flow in the m direction; іп? /ѕес. 
journal radius 


radial and tangential reference axis 


fluid velocities in the x,y, and z directions res- 
pectively 


perturbation of journal center; in/sec. 


film width at any point measured from the circunfer- 
ential groove; in. 


horizontal and vertical reference axis 


film coordinate axis 

а 

dt 

width of a oil striation at the downstream boundary 
of the cavitated zone; in. 


width of a oil striation at the upstream boundary of 
the cavitated gone; in, 


angle between bearing load and line of centers; ra- 
dians 


rotational speed of journal; radians/sec. 


relative tangential velocity between journal and 
bearing, (W- 28); radians/sec. 


AL La 


dimensionless squeeze, 
Pa © 


viscosity; lb-sec/in?. 


angle between line of centers and any point on the 
journal; radians 


Shear stress; psia. 


vi 





INTRODUCTION 


Background 
The knowledge of the dynamic properties of a journal bear- 


ing is necessary in the design stage both to ascertain the 
inherent stability of the bearing itself and to predict the 
overall response of a bearing supported system to external and 
internal excitations. Some of the many problems which can be 
encountered are self excited half speed whirl, excessive am- 
plitude of vibration of the journal within the bearing near 
the critical speed of the system, and the inability of the 
system to recover from a dynamic load. 

Since Reynolds‘derived the three dimensional, continuity, 
differential equation for determining the pressure distribu- 
tion in an oil film between two moving surfaces, all journal 
bearing investigations have used Reynolds! equation, usually 
in one of two simplified forms because of its mathematical 
formidability. In the two dimensional or infinite bearing sim- 
plification, the pressure induced flow in the axial direction 
is neglected. In the short bearing simplification, the pres- 
sure induced flow in the circumferential direction is neglected 
by assuming that it is much less than the velocity induced flow 
caused by the relative motion between the journal and the bear- 
PAS. 

However, defining and applying the proper boundary con- 
ditions has proved to be the most difficult task in obtaining 
a good solution to Reynolds! equation, even in one of the ге- 
duced forms and under steady state conditions. Sommerfeld, in 
1904, was the first to obtain an exact solution to Reynolds! 
two dimensional equation, but his solution allowed negative 
pressures in the diverging portion of the bearing. This led 
to the theory that under certain conditions of loading, vapor- 
ization could occur because a fluid is unable to support any 
appreciable tensile force. This theory has since been verified 





by making observations of tranlucent bearings which clearly 
show film rupture(1],.[2].[3].[^]- Theoretically, the cavita- 
tion pressure should be the vapor pressure of the fluid, but 
in actual practice the cavitation pressure has been found to 
be very nearly ambient pressure because of entrapped air and 
other gases. 

In 1952, Ockvirk[5] developed a steady state solution of 
Reynolds! equation using the short bearing approximation and 
assuming that the entire unloaded half of the bearing was ca- 
vitated and at zero pressure. More recently, Fedor[6] and 
Donaláson[ 7] obtained solutions of Reynolds! three dimensional 
equation using series solutions and applying Sommerfeld's 
boundary conditions, but requiring that the feed pressure be 
sufficiently high so as to suppress cavitation. However, this 
is an unrealistic situation in most bearing applications, par- 
ticularly those in which stability is a problem, 

Two recent analyses of unsteady motion are of special 
interest in this investigation. Holmes[8] extended Ockvirk's 
solution to the particular case of a rigid shaft supported by 
symmetrical bearings. Richardson[9] also used the short bear- 
ing approximation, but defined the cavitation boundary more 
rigorously by assuming that there was & zero pressure gradient 
along the boundary and by assuming that the shape and extent of 
the cavity was a function of feed pressure and squeeze..Rich- 
ardson obtains a solution to his analysis by digital methods. 
However, both Holmes and Richardson violate continuity since 
they do not consider the oil entering and subsequently leaving 
the cavitation zone. 

In 1957, theories which took into account continuity of 
flow through the cavitation zone were independently developed 
by Floberg[10] and Rightmire[11]. Experimental work by Flo- 
berg and Ishii[#] has indicated a close agreement with the 
continuity theory under steady state conditions. 





Purpose 
The purpose of this investigation is to analyze the un- 


steady motion of a circular, circumferentially grooved journal 
bearing, using the short bearing approximation and the bound- 
ary conditions of continuity. <A circumferentially grooved 
journal bearing was chosen because its line source makes the 
feeding arrangement independent of the journal attitude within 
the bearing. The results of this analysis will be presented 
in two forms: a computer program which can be used as a sub- 
routine in a systems analysis, and in the form of tabulated 
data.of spring and damping coeffecients for use in hand calcu- 
lations. 


ANALYSIS OF UNSTEADY MOTION 


Assumntions 
Invarrivings ас а solution of Reynolds' equation which 
Will adequately describe the unsteady motion and dynamic pro- 


Dervies of a circumferentially grooved journal bearing. the 
Following assumptions will be made: 


tine fluid is Newtonian and incompressible. 

EE ne viscosity is constant. 

3. The flow is laminar. 

4. Tre mass of the fluid is negligible. 

5. Surface tension is negligible. 

6. The fluid adheres to the journal and bearing surfaces. 

END ere 15 пс distortion of the journal or bearing 
O O nyaccodynento pressure on iy lehnt 
and moderately loaded bearings will be considered, 

Freche fluid film iS so thin that it has no curvature. 

СУС ус ОО ТОЛ йаг ореге пе NO pressure 
pradient across (һе туш. 

ШОО The vapor pressure of the fluid is equal to ambient 
pressure. 

11. The pressure induced flow in the circumferential di- 
rection is negligible compared to the velocity 
Dnguceserlowsrn Тһе crrcuaterential direction 

12. In the ceviteted region, the fluid divides into an 
infinite number of streamers. 

The last two assumptions will be discussed in greater 
detail in later sections; the remaining assumptions are self 
explanatory. 


ла! Bearing Oil Film Theory 
Figure I illustrates the general case of the unsteady 
Bon ol a journal within a fixed rigid bearing. At any 


Gnsvant in time the journal has a rotational speed W about 





шіге I 
Journail-bearing motion referred to a fixed point in space 
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Figure Il 
Eournal-bearing motion referred to the line of centers 
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Pigure ү 
Oil film forces in the x direction 
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Figure IV 
011 film forces in the z direction 
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meso Onn center CJ, an orbital speed of Й about the bearing 
center CB aná a radial velocity e. However in analyzing tnis 
Byevem, it is more convenient to fix tne line of centers in 
Ecce end use that line es a reference as shown in Figure II. 
Then the journal has a clockwise rotational speed of (W - Y) 
eng the bearing has counterclockwise rotational speed of й. 

Kc forces acting in the x ana z directions on an ele- 
metal cube dx dy dz are shown in Figures III and IV, where 
uU, v, апд и аге the fluid velocities in the x, y, and z di- 
wee ons respectively at any point. Equilibrium of forces 
Оле x and y direction yields: 





до ER 5С, P 

òx ay. sz UL) 
ор | 9U., Г: (2) 
¿Z Ay ax 


meen tne assumption that the film is so thin that has no 
pressure gradient across it; 


др. 
зу 509 


The sheer stresses are obtained from Newton's law for laminar 
viscous flow[12]: 


Substituting these expressions into equations 1 and 2 gives: 











др 2 А2 1 + ae A 

ELE RE иа [ды аы) 

à x dy ӨУд% ““ӛуі ду д2 | (3) 
& D ° = Y 14 2 ur : È 2 W а W | 

— = 23 — = Af —= — Im no (4) 
фа ay az A Sr tor 


again considering the fact that tre film Is Very Cade an. the 


Ж ЕЕ ТЕсСсіоп, so that the velocity gradients in the y direc- 


tion are much greater than the velocity gradients in the x 


and z directions: 
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equations 3 and become: 


Be E 
EB ^ 
A 
за И 


S (5) 


оғ 
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L 
Фи 


ду 


(б) 


Withthe assumption of constant viscosity, equations 5 and 6 


can be 
_ 1 
E 22 
E 
W = 2 


integrated twice: 


AD 

m y? + Су + С. (7) 
òP 

SZ у + Суу + С (8) 





Neo 52255, алп Os are constants of integration. 
From Figure II, the boundary conditions are: 

Sey = Of: = ш @)В, и = 0 

at y = h: u= - #В, и = O 


where h is the film thickness at any x. 


Applying these boundary conditions to equations 7 and 8: 





NE Ox Eo uu ү 
| 1 àP 2 


With the assumption that the pressure induced flow in the 
EEUcumferential direction is much less than the velocity in- 
duced flow due to the relative motion between the journal 


and the bearing, equation 9 becomes: 


uz - 28 + (WW. G)R ГІ) 


The boundary velocities in the y direction are: 


1 
O 
«1 
1 
e e 


at y. 
s пура ие 


IE continuity: 


4 
| 
е 
| 
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Euostcrtuting equations 10 and 11 into the continuity equation: 
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Integrating the above expression from y = O to y = h and 
applying the boundary velocities in the y direction, the 
final result is the short bearing approximation of Reynolds! 
equation for unsteady motion: 





NE S m, ша 
e > ns d (12) 


where A = (W)- 2) 


This equation is much easier to manipulate and solve 
than Reynolds'! complete equation, and for moderately loaded 
bearings with length to diameter ratios of less than 3 the 
feos Of accuracy introduced by tne short bearing simplifi- 
Cation can be expected to be of the same order of magnitude 
as some of the other assumptions which were necessary to 
make. A comparison between an exact solution obtained by 
Donaldson [7] under steady state conditions (6 = 0, ё = 0) 
and the results of this investigation is presented in 
Appendix A. it is proposed that the same deviations could 
be expected in the case of unsteady motion between exact 
and approximate solutions. 

Again assuming constant viscosity, equation 12 can be 
EN rated twice with respect to Z to obtain the pressure 
distribution on the bearing surface: 


10 





pesen S (CC On tanto mor integration to be dez 
mermined by сле film boundary conditions. 


A pressions Tor the nnii flows'in the x and z directions 
will be required in the next section for determining the cav- 
tation boundaries. These expressions can readily be obtained 


by integrating the velocities u and v across the film as 


FOLLOWS: 
a 
Q, = u dy 
° 
n 
Q и у ду 
Jo 





уз ёр | 
dx == 127 dx E 
е. Bh и? sp 
8- - (Ш- 26) 2 7124 52 (TS) 


Film Boundary Conditions 


The unwrapped oil film for one half of a circumferen- 
tially grooved journal bearing is shown in Figure V. The 
approximate shape of the cavitation Zore as observed experi- 
mentally [1], [2]. [3]. [^]. is defined by the boundary 
l-2- 3; for steady motion, point l occurs immediately 
downstream from the point of minimum film thickness (x= mH) 
where the rate of film expansion is the greatest, but in the 
case of unsteady motion it must be assumed that cavitation 
could commence at any x. Since the vapor pressure of the 
Tluid is assumed to be ambient, then the entire cavity must 
be at ambient pressure, with the result that one of the 


n 





boundary conditions is: 


aç Z = W, p = Da 


where w is the film width at any x and is egual to L 
where cavitation is not present. 


Boplying this boundary condition to equation 13, tne result is: 


2 


——— — - —y(z2° - w?) +C (z =- w) (16) 


Pigure V 
The unwrapped oil film 
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With the assumption that the fluid adheres to the sur- 
faces of the journal and bearing, the fluid must split into 
one or more striations upon reaching the boundary 1 - 2 in 
order to satisfy continuity and to prevent subambient pres- 


sures from occurring. This is shown in Figure VI for the 
particular case where h,>n,. 
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Figure VI 
Continuity of flow across the cavitation boundary 
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The upstream boundary 1 - 2 is determined by considering 
tne flow througn an elemental volume h dx dz at A. Then the 
flow entering the boundary at A is 


9, - Q,dx + 9,02 


IN 
From equations ll! and 15 this becomes: 


h? aD 


кё _ 929 A 
Qu, = &(U- 22)Rhdx-,27 ұс 42 


Since there are no pressure gradients in the cavitated zone, 
the only flow is velocity induced flow caused by the rela- 
vive motion between journal and bearing. Thus the flow 
leaving the boundary at A must be: 


дш. = E(0 - 2P)Rh dx 


Therefore, by continuity of flow across the boundary: 


13 
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along boundary l - 2. 


Applying this boundary condition to equation 16 to evaluate 
the constant C, the pressure distribution in the region 
bounded by z = 0 and 3-4 - 1 - 2 is found to be: 

h 
= 2L tg A98 
A = Ox 


~ 


A - =) (17) 


where p' = Т. - Pa 


The first term on the right side of equation 17 is the 
hydrodynamic pressure distribution and the second term is the 
hydrostatic pressure distribution. The slopes of these two 
terms at 2 = н, are, respectively: 





| wAf oh : 
ta = (RA - 2y) (18) 
0% hd Hn" ax 
9р р! 
22 hs E на 


tne existence of cavitation at any x upstream from point 2 and 
downstream from point 3 can be ascertained by letting w = L in 


the above expressions. Then the condition for cavitation to 
eur is: 


я ' 
ar -2y) > (E=) (20) 





15 








i.e. the positive slope of the hydrodynamic pressure is 
greater than the negative slope of the hydrostatic pressure. 
Mais Ts illustrated in the left diagram of Figure VII. 


Figure VII 
Axial pressure distribution in region 1 - 2 
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equation 20 (unstable) equation 22 (stable) 


If the condition for cavitation is satisfied(equation 20) then 
w adjusts so that: 


Ж, (El 
СЫТ 7 |22 hs 


Substituting equations 18 and 19 into the above expression, 
tne film width along boundary 1 - 2 is 


prune 


3M (RA -2y) 


2 


W (21) 
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Making tne substitutions: 


n сп 
95 = = = sine 


Уш - é cose 


Equation 21 becomes: 





= p'h? 
E плав 
34 (enAsind- 2& cos®d) (22) 


Hf equation 20 is not satisfied, then: 
pus IL 


It can be seen that equation 22 is invalid for the down- 
stream boundary 2 - 3 by considering the steady state case. 
Then equation 22 becomes: 


=p he 
“= [34 enWsine 

Under these conditions, the expression under the square 
root becomes negative at angles greater than 8 = 21, implying 
that the cavitation zone cannot extend past 9 = 2m. Physically 
this would not necessarily be true for high journal speeds and 
Mew teed pressures. 

The downstream boundary is determined by considering con- 
tinuity of flow through the cavity as illustrated in Figure.. 
VIII. Let the width of one striation at the upstream bound- 
ary be Az. Then the amount of fluid in the striation as was 
previously determined is: 


3ARBh* Az 


16 





where h* is the fílm thickness at z = w on the upstream 
boundary. 


Figure VIII 
Continuity of flow across the downstream boundary 





With the assumption that the fluid adheres to the journ- 
al and bearing surfaces, tnen by continuity the same amount of 
fluid must leave the cavity at Zz = w at point B on the down- 
stream boundary. Since we have assumed that there are no 
pressure gradients in the cavity, there can be no motion of 


the striation in the z direction. 
Summing up the total flow in the unit volume h Ax Az 
at B, the result is : 


h? n 











SARNA z - 12H 3z vx | Ax = 3ARh Az (24) 
ia 
From ecuation 17: 
o = _ mes (25) 
ü а; 
ұу ¿ 


ши 





Substituting equation 25 into 24: 





&ARn* Az -= pu (oh ES ES P Ax = $ARh Az 
124  n? ох W 


Assuming an infinite number of striations, and then taking the 
limit as Ax—dx and Az—dz: 





t ) 
=ARh* dz | ERAS ey) = d dx= i.ARhdz (26) 


Making the substitutions: 


ах = ваө 
E! = dW 
iw 
oh. - 2 sine 
y s - e cos 0 


Equation 26 becomes: 





3A Rh*dw m LE enfisin8e+ 2e cos 8) в de 
J 
= ФЛАВИ ам 
or: Пере 
dw _ ¿w(enÁsin 6 - 26 cos 0) - ÓAf w (27) 
E A(h -= hë) 


The film width, w, in equation 27 is in terms of two vari- 
ables, dw/d® and h*, and so must be solved by an iteration 
method. 


18 





013 Film Forces 





K О i MeTorcesrwnicn will be considred are: 
l. Hydrodynamic, hydrostatic and shear in the complete 
film région pounded by z= O and 4 - 1-2 2 3 = 1 
in Figure V. 
2. Shear in the cavitation zone. 
The differential hydrodynamic and hydrostatic force on a 
 ГІр сі the journal w, В 4015: 


ги 
aF, ра аө 


220 
 егітигілг іп p from equation 17: 
hi һ > г. | 
т? Е ЕЛЕЕ паа 2 _ тут ДЕ m 
ar, = || à (RÍL = туг Zu) + p'(1 =) | az R ao 
{ 


55 «d 
мез | 


wnich when integrated becomes: 


y Го и в Ah T 1 
=! - - 


M | 2n: ox 





Resolving this differential force in the l and 2 directions as 


shown in Figure IX: 


E. = àF, сос 0 


d Fpa = ағ; SM 6 


“ne tangential shear stress in the complete film region 
due to the relative sliding velocity is (see Figure ITI): 


2 е е 
T= wRr- fe cos O9 + esino) 
However, Be< WR anà é& (DR, so that: 


. Ма) 
Gx R 
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A En 7 


E = 
| f i 
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Fr d 
ы 





Figure IX 
Coordi Wave cyscven for oil film forces 





Then tne differential force due to shear on a sirip и, Rdeis: 


dF. = C wR åo 


5 
Substituting in the value for D: 


айын, 
© h 


at 
Resolving the shear force in the 1 and 2 directions: 
ак, = = aK sinodo 


Alea = dE. cos с 96 


In the cavitated zone, there is no hydrostatic pressure, 
Since p=p,. Also, since there are no pressure gradients, 


tnere is no pressure flow and thus no hydrodynamic pressure 
exists. 
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The shear forces in the cavitated zone are obtained by 
assuming an infinite number of striations as in the preceding 
section. Refer to Figure X where one striation is shown. 


Figure X 
Continuity ог Single oil striation 
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con tinuity: 


F 


2 ARh*$* 2 LABh6 


е 


where &* and $ are the striation widths at the upstream 
and downstream boundaries respectively. 


Then: 
h*$* = he 
or: 
9% 
>= h 


oA. 





Substituting in the expression for $ from above: 


а 
АШВ п 


at. = 
i: h n 


In the limit, as $* — dz: 


UR? 


af = 


h*dz de 


Then the differential force due to shear across the complete 
width of the cavitated zone at any O is: 





а 
Af (UR? 
аж. = WR? ne az ae 
2 ) h 
Š = vV 


where the integral must be evaluated numerically. 


Resolving the cavitated zone shear force in the l and 2 direc- 
ESOS: 


due == d sinode 
әс 


о. 
txj 

e 
|] 


ар. cos 8 d9 


cumming up the forces in the 1 and 2 directions and 
EUDpStitutions: 


ЕЗ 

(t5 

ж 
[e^ 
b 
GI 

cet 
es 
Gp 


Then the differential forces in the 1 and 2 directions are: 








з 
W R В t 
EE 4 (Acnsin9 - 2ecos @)cos @ + zn cos O 


217 И 2 ) 
: ne 
_ MWR W sine MAR Т 


- 
* " sin 6 | de (28) 


2 





| wR - 8р! 
dF, =| n onmsinessc2ecosO9)sin6- 
| 2h 2 





sino 


ARÈN Mivr*I ] 
los uoo (29) 


where I - | h”dz 
ұу 


F =| åF 
Jo 


whicn must be evaluated numerically. 


The attitude angle is 


2 
g = Бат (ро) (30) 


Spring and Damping Coeffecients 


Refer to Figure XI. At a given instant of time the jour- 
nal center is at CJ and is being acted upon by the forces F, 
and F, which were evaluated in the previous section. Now the 
journal is given a perturbed displacement, d, and a perturbed 
velocity, v, to the position CJ'; then the forces acting upon 
the journal are of the forn: 


К ткыш с пак в те) 


T 
1 


io 


Fa 


m 


pou) pus) 


where F,, and F,,are the steady state forces at CJ; and K 


and B are spring and damping coefficients. 
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Figure XI 
Perturbation of journal center 
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In evaluating the previous equations it is first nec- 
Eeeary to transfer the forces P, and F5 to the r and s axis. 
However, it is more convenient to continue using differential 


Boroces; thus: 


ar 


I 


^ AF cosx = dF sina 
dF, = dF,sinA+ dF cos œ 

Since a small perturbation is being considered: 
eo CJI << CB CJ! 

men: 

S 


nn S 
S e o — р a 
in AZA= e = C 


cos A al 
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With these approximations, the differential radial and tan- 


gential forces become: 


nw 3 (31) 
r | z cn 
S 
dF; = dF — + dF, (32) 


substituting tne eauations 28 and 29 into equations 31 and 32: 





Ми в ° | 
ЧЕ, = E cos 6 | - 26 cos в + Ane sine|+2 ; piwR cos Ө 





2h 
ғ” 2 (21 A —- 2 — ; 
A = W sine — + sine s (en B cine 


+ zpwRsin® 


Па 22 


a : 
|- 2e cos 6 + A ne sino 





- ^ 
ы IR? А Hj ck 
Ыы ы — T cos ө\$ ае 
h y? / 
(в S ЖҮЗЕ В e ° 
аР. = | са аа соз Ө| - 2е соз 6 +. сп зіп Ө | + 2 риг соз 
iL) Re MUD RWI E yy 
BAe Wen ы а E S sin © 
h ny 2n? 
- 26 008 9 + Ano sino) + 2pwhRsine 
MU RE: Aun? T | 
RAN у HUS eS ef де 
n y? 
Let: 
E Cn 
h = с(1 + псоѕ ө) 
ПЕЕ cs 


cn 


Amos tituting the above into the equations for dr, and ағ, 
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and rearranging: 





: CA м Я Г 2r Cose W ne sino 
dF. = --- сов Ө (- mm t 
r 205 Ве сова ПЛ сос С) 
2ssin6 |, j cose . Row sine 
(1 +ncos 8)? | “ P (1 +n cos 6)c 
ДД ВРТ sine 4 w?Rsine Г 275 
e2(1 +m cos 8)* - 263 | сп(1 + псоѕ Ө)? 
(iis sins 255 віп Ө | p wHs sin © 
(1+ncos 0)* ~ en(l+ncos 86)? | Den 
AU)BR?ws cose  ALuR?IscosO > 
лы е еи (33) 
c"^n(l-ncosG) c?n(1-4ncos off 
Е M wR cos Ө 2rs cos @ ñ lis sine 
s 2c? (1+ псоѕ Ө)?есп (1 +m cos 6)° 
255 біп Ө | p'swBeosOo  J4pR?ws sim 6 
сп(1 + п соз Ө)? | Zen c2n(1 +n cos 6) 
AU R?Is sim € ÄH wR sine f 2г соз Ө 
c?n(1-ncos 0)* 26 L (1+ncos 6)* 


Uy nc sin e 2s sin O 
+ > рив зи е 


(1 + и соз 6)” j (1(+псо Ө)? 


Ala?weos6 MiUR?Ic х 
E NE сос О E (зц) 
с(1 + псоѕ Ө с2 (1 + п соз Ө)? 


- 


Equations 33 and 34 are non linear inn. However, with the 


assumption that the perturbation is small, the equations can 


be linearized. From Figure XI: 


cn, -r - cnooscícon 


Or: 
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where E << папа 15 time dependent. 
By Taylor's theorem: 

fín +) = f(n,) + 2£'(n,) 

neglecting second order and higher powers of 5 5 
Then, by also neglecting other second order terms: 


S UE с 
E Гүл т INS fin) 


D fín, +2) <3 гҒ(п,) 


T I: mL 
с т (юй, br) p (л, ) 


О | 


2 


r 
S ac COS fino 


rs= 0 


ІД 
y 
! 
O 


Linearizing the terms of equations 33 and 34 as above: 


2r cos Ө E 4 en sine 
(1+ псоѕ Ө)? (1 +псоѕ Өө)? 


(Шт ѕіп Ө (1+ псоѕ 6 – Зп соѕ Ө 28 sin 6 | 
КА — RÀ XA - a 
1 + д соз 6)3 (1 +m cos Ө) 
Alun?w sino  , MUR?wrsinO cose 
RM 1M EE ouo o orm Me y sa 
с(1 + й сов @) с? (1 + псоз 0)3 


AMUJR?ISsine А 2MUIR?Irn sin 6 cos 6 
с2 (1+ псоѕ ө)? сэ (1 + псозѕ e)? 


| Ai WÊR 
ao, = соз Ө | - =- 





r 263 


+ Е ривВ соз 0 - 
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X wRsWsin®? 6 pwRssíinO AU R?wscose 


20° (1 + п соѕ ө)? zen E сп О + п соз Ө) 


HW R?Is cos O 55 
cén ( +m cos 6)2 


(35) 


Др. = Ж w?RsUcos o sin os p swR cos 8 AUI ws sino 
> 2c? (14nocos 8)? 2cn c2n ( 1 +m cos 0) 


М Ц) R2Ts sim 8 | Пит Н зп Ө d 2r cos © 
с?п (1 + псоз 9? 2c? m (1 +m cos 8)? 
я Unc sine UE sin6(1+ncos6 - 3ncos 6) 
(1 +n cos 9)? (1 +m cos 6)“ 
25 ѕіп Ө A Lidl) Rew cos 6 
er zn боол о 
(l+ncos6)* | c(1+ncos Ө) 
MUUR? wr cos? 6 MU R*Icos © 
с? (1 + псос 6) 2 с2 (1 + псоз 0)? 
2 4UR?Ir cos? ө) 
- — h --. (36) 
с (1 + йсоз 0) y 


Equations 35 and 36 are of the form: 


CP. = - dK,,r - dB,,r - dK,.s + d8,¿S + di, 


where: 
МОм? В 
КИЕК 222 eee COS OL neos O- Зисоѕе 
rr E ETT 3 ) 
ALUR? w sin Ө соз Ө 2ЛМА)В21 51п Ө соз © а (37) 
ce?(1+ncos0)? "^  c?(1-4ncosO)? 
prs 2 
EE кее Ө (38) 
Ec cos e) 
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ДШ Ват соз е 


P'WR sin © Ñ 
e2m (1 + mn cos Ө) 


ак „АЕ sin” © 
rs 7) 262 (1 + псоз 0)” en 
АЛ 2 5 
a I cos @ MT (39) 
с?и (1+ псоѕ Ө)?) 
(40) 


yy 
ав, = em к ае 
r с” (1 + а соз е) 
PWR cos е, MU R?w sine 


A4 w? RU) cos 0 sán 6 
са (1 + псоз 6) 


К = 
$° 2c? (14 nocos 0)? 2cn 
Adil R?I sine 
- de 41 
c?n (14 ncos 6)? S 
ÆA ни В sin? 6 
ав... = Е — 40 (42) 
эз c*(1+necos 8)’ 
ак... = АНИ нов sin? O (l+ncos6 - 3ncos 9) 
sr 204 (1+nc0s 0) 
{uU n? ° 2HWR*T cos? © 
sA R? w cos? e ES Se. dio (43) 
с? (1 + исоз 6)” с” (1 + псозв Ө) 
Ма В sim бсоз 0 
ав = - -- де цц, 
sr e*(1l +n cos 8)? n 
then: 
2# |7 
L. | ОК num 
а 
21 
Ba ж E dB nn 
© 
wnich must be evaluated numerically. 
The final result is: | 
3 + F (45) 


mr p Y K съв SS 
Py rr rs rs 


Ша = 
Д 
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а Е 2 ° e 1 
i K.,S - BS tK rT +B Ur +F (46) 


However, it is generally more convenient to work with 
Ú tesian rather than polar coordinates. Therefore the in- 
tegrated forces, and spring and damping coefficients will be 
transferred to the x - y axis as shown in Figure XII. 


Figure XII 
Journal motion in x - y coordinates 
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From Figure XII: E 
r=ysinf +xc0s f (47) 
s = ycosf -xsinf (48) 
E y cing +x cos # (49) 
s-ycosf -xsinf (50) 
Е =F,Cos Y - F¿sinf£ (51) 
Б = Р. ѕіп + Е, соѕ 0 (52) 
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Equations 47 through 52 are substituted into equations 45 and 
46; 


Ше = - ЕУ sing -*xcosf)- Bay sin Ë + x cos f з 


- Ke (ycosf -xsinf)-t Bu (ycosf -xsing) (53) 
Ен - Kgs (y 0088 - xsinf )- Bes (y cos Ø - x sin f ГЕЛ, 
Их сову) = Be, (ysinf + xcosø ) (54) 
Substituting equations 53 and 54 into equations 51 and 52: 
F = - К. sing cosÜDy - 5. 45652 7 X - В,, віп 2 cos йу 
- EOS @ x Ke cos? Y y + K re sing cos fx 
cos? 8 y - B 


+ Bye 


ee sin Y cos f x ЕС Е05 0 
+ Kos cos sinf y - Kee sin? f x + B sg cos ø sin y 


- Вс; sin? f = K., sin? f y - Ks. sin f cos £ x 


- B, sin” f y - Bon sin? cosßx - Pe sing 


6^ 


Fu = - K,, sin? f y - K, sing cosgx - B,, Sin? f y 


p Sing cos f x E са Й cos dy + и sin? f x 
+ BL. sing cos fy - В... sin? f x + к сат. 

- K,. cos? f y + Kss sing cos ú x" E cos? f y 
EB. sin cos g x ЕЕ біп соѕ Йу + Ке, cos? 0 x 


Ber sin Y созбу + Be, cos? ff x + a cos Ø 


EL 





Then: 
P. (Хх 25-4 K ¿y Y - Buy Y en 
Fy = - Ку,у - By,y + Кух - В„„х 
Where: 


Ku = Kp, cos? Ø + (K..- Be) сз сос ПК се 


yy = Kss cos? Ø + (Kp¿- Kgp)sinfS cosf+ Kpp sin? Y 


PS 
x 
| 


y= K rs cos? f + (КК. - Ку) дай соз + К., sin? f 
K,¿= Kgp cos” + (Kgs- Kpp)sinf cos É + Ky. sin? f 
Bix = Bpp cos? + (Bho + Bep) sin@ cos $+ Beg 5117 
; cos? f - (Best Bor) sing cosp + Bpr sin? Ø 
Byy = Bep Sin? Ø + (Bpp- Bg¿) cos f sinf- Bys cos” Y 


Byy = Byy 


Computational Procedure 


(56) 


(57) 


(58) 


(60) 


(61) 


(62) 


(63) 


(64) 


The general case of unsteady motion of a circular cir- 


cumferentially grooved journal bearing has now been analyzed 

and the differential forces and differential spring and damp- 
ing coefficients are in forms wnich can be solved numerically. 
A detailed procedure, using a high speed computer, is given in 


Appendix B. However, the general computational procedure is as 
follows. For a given bearing geometry (І,, В, С, апалу) апа 


given operating conditions (Ww, Й, n and e): 


Be 





1. Determine direction of flow (+ Ө) and locate up- 
stream and downstream cavitation boundaries. 

2. For each ©, calculate the film width. 

З. For each Ө, determine the differential forces and 
differential spring and damping constants. 

4, Integrate numerically to obtain the radial and tan- 
gential forces and spring and damping coefficients. 

5. Calculate the bearing load, load angle, and horizontal 
and vertical spring and damping coefficients. 


Amplified Analysis 
It should be noted that equations 37 through 44, and, 


consequently equations 57 through 64 are independent of 2,0, 
de, aná AS due to the fact that (see Figure XI): 


_ lim 0 


Б = do dd (ак, 4) 
Cy 
dB „ e EE 


> О àv “р^ 
e 
Since dK,,, and dB,,, are functions of w, then w must also be 
independent of й and è in the above mentioned equations. 

If now frictional effects are neglected on the assump- 
tion that they are small compared to squeeze effects, then ~L 
and U) can be used interchangeably througnout the entire anal- 
Sis by virtue of tne choice of reference axis. Then the equa- 
tions for the spring and damping coefficients and the equations 
for the steady state forces Го апа Ir can be non-dimension- 
alized and expressed in terms of just three independent vari- 
ables; eccentricity ratio, n, pressure ratio, PR, and dimen- 
sionless squeeze, /X. Let: 


ry 
N 
ojo 
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Then the vertinent dimensionless equations are: 




















ABs @ 


p LR 


(PR - 1) (14 nocos 9)? (65) 
Se ca oS т 
(Ey 2 1 КЕССЕ С) (ЕЕ L] (66) 
L' 2(cos 09- cos €*) ónA*(cos 6 - cos 6) 
sínOcos O 
EM at E *(*) (PR- 1)eose de (67) 
y n? © 2 
NO Eum E UE + (= -) (PR- 1)sine àe (68) 
J® wo. sin® cosé 
= = I Genese (1 - ncos O0 - За СО5 6)ds 
(Sy) 
=| 2 S “а уй Асо56) г 2n e = ) г) sin 4 (70) 
| A сіп Өсов Ө 1 үт I 
=:— — (y 4 (PR = 1 — 0:40 
| 2 2° (1 + ncos 0) ТА D es 1 
(72) 
Ле sin” © | 
E S ix conu + uncos O - a соз ©) (72) 
NW. cos? © a (73) 
(1 + ncos O9)? Б 
I sin? 6 
= Дж de (74) 


(1 + ncos 0)? 





AB,r C sino cosO 
o LR (1 + ncos e)” 


ae Z 


The general computational orocedure is the same as before, 


Серт that all frictional effects are neglected. 
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RESULTS 


Figures ХІІІ, ХІҮ, ХҮ, ХУІ and XVII illustrate tne shape 
and extent of the cavitation zone for various combinations of 
journal motion obtained by this investigation. Results obtained 
by Ishii and Richardson ( the latter in dashed lines) are also 
presented Yor comparison. The direction of the circumferential 
flow in Figure XVI is the reverse of the circumferential flow 
ШІ Гірпге ХУ, thus Figure XVI is a mirror image of Figure XV. 

Tables l and 2 show a breakdown of the effects of circun- 
Terential squeeze, radial squeeze, feed pressure, complete film 
friction, and cavitation zone friction on the forces and spring 
and Gdampins coefficients for the examples in Figures XV and XIV, 

Levulated data for the simplified analysis which negiected 
Srictional eficts is given in Appendix C for numerous combina- 
pons Of dimensionless squeeze, pressure ratio, and eccentrici- 
ty ratio. The parameters were chosen so as to cover as wide a 
гапве of bearing epplicavions as possible. 

Table 3 is & comparison between the data in Appendix C 
end data obtained by Holmes E who analyzed the vibration of 
a rigid shaft by extending Ockvirk's analysis and boundary con- 
ditions. The data in Appendix C was converted to Holmes’! nomen- 
clature as follows: 





ho 
зо 
E 
B 
! 
= 
E E 





Wes = вщ ан 
por nm EI 
Wnere: 
^ ^ | 3. 
HE Pome den ee ton JE 





marca 1s a Tora ot Оскузги с loed number. The factor of 
2 із necessary beceuse tnís investigation uses a circun- 
ferentially grooved bearing with a total length of 2L. 


Since Ockvirk's theory defines the attitude angle for a given 
eccentricity ratio, it was necessary to fix both 9 and n in 
fae comparison. 

The program in Appendix: B is written in a manner so that 
ИШ/ ог ell cf the values can be called without any major changes. 
Frictional effects can be neclectec simply by deleting the in- 
dicated cards, thus reducing the program to three dimensionless 
variables as in Appenix C. It should be noted that for sim- 
ра: Бу оп! у one half of the bearing was considered in the 
analysis, but that the program sand tabulated аба are in terns 
of the entire bearing length 2L. 
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BESCUSSION OF RESULTS 


For the steady state exemple of Figure XIII, Isnii's results 
иам ест но озе to the resuits olwthis investigation, consider- 
ШЕ лас while Ishii used continuity of flow through the cavity, 
he averaged the flow leaving the boundary 1 - 2 end then used 
this average to determine the boundary 2 - 3. Thus it should 
be expected that Ishiits point C should be outside the bound- 
ary of this investigation since the average flow at Ha is less 
than tne straight line flow, while point D should be inside the 
 пайату of this investigation since the avearage flow at Wy 
is greater than tne straight line flow. However, h does not 
Vary appreciably between point 1 and 2 due to the large slope 
of the upstream boundary, anā thus the average flow through- 
out the cavity is very close to the exact flow. Also, some 
minor disagreement could be expected between machine and hand 
Baleulations. 

The deviation between this theory and Richardson's is 
quite pronounced along the downstream boundary, particularly 
in the steady state cases of Figure XIII and AIV; closer agree- 
nent is reached as tne ratio of circumferential flow to axial 
flow is reduced as in Figure XV and XVI until there is perfect 
agreement іп Figure XVII where only axial flow due to radial 
velocity e is present. None of the sample runs indicated any 
differences along the upstream boundary; tnis was to be ex- 
pected since both theories use the same method for calculating 
the upstream boundary. The values for load and attitude angle 
are surprisingly close in Figure XVII, considering that dif- 
ferent numerical integration techniques were employed for the 
two theories. Thus any deviations in the results can be attri- 
buted directiy to the differences between the downstream cavi- 
tation boundaries. This is best illustrated in the steady state 
examples of Figures XIII and XIV, where the most pronounced 
difference is in the fourth quadrantí"(1.57<0<2m).. 
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Richardson's theory yields a greater film extent in this areca 
and consequently a greater film force in the direction of the 
applied load, thus reducing the net bearing load and increasing the 
Load angle. The samo analysis can be applied to the other 
examples although it is somewhat more complicated in the dyna- 
mic cases, Richardson did not compute the spring and damping 
coefficients and so no direct comparison of these can be made, 
but by comparing the bearing loads and load angles for the 
,Gynamlo cases,, it can be deduced that the same order of mag- 
nitude of differences in the spring and damping coefficients 
would exist by virtue of equations 45 and 46. 

The purpose of Tables 1 and 2 is to show the small effect 
that friction has on the film forces and spring and damping co- 
efficients, and therefore to justify neglecting friction in the 
preparation of the dimensionless tables in Appendix C. In gen- 
eral, friction can be neglected in all cases except when 
(LU - 20) and e are very small, as in half speed whirl. 

Table 3 illustrates the effect of the choice of cavitation 
boundary conditions on the spring and damping coefficients. 

Using Ockvirk's boundary conditions, Holmes assumes that the 

entire half of the bearing is cavitated , an assumption which 
approaches this theory only at the higher eccentricity ratios. 

Thus there is close agreement at n = .6, while at n = .2 the 

two theories yield somewhat different values. At first, it 

would be expected that this theory should always give less 

values of Ese Kos’ Bor’ and Bos at low eccentricity ratios 
since these coefficients go to zero in the completely uncavi- 
tated situation(see Appendix C), but this is not so because 
the film and the film forces obtained by this theory are un- 
Symmetrical about O = 1.5 , and thus a clear comparison can- 
not be made. However, by comparing the results of this inves- 
tigation for a high and a low feed pressure, it can be seen 
that as the feed pressure increases causing the film extent to 
increase, the radial spring coefficient decreases. 
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Sternlicht| 13] used Sommerfeld's boundary conditions and 
digital methods to solve Reynolds! three dimensional equation 
end to obtain spring and damping coefficients for a full jour- 
nal bearing. However, no quantitative basis for comparison 
could be found, although it is significant to note that Stern- 
licnt found the damping cross terns Bxy and Byx negligible, 
while this theory and Holmes! analysis found these coefficients 
to be of the same order of magnitude as B and В." 

The general validity of the results can best be ascertained 
by examining the initial assumptions that were made. Since 
low to moderately loaded bearings are considered in this investi- 
gation, bearing distortion and viscosity variations due to pres- 
Sure can be neglected. Also the choice of a circumferential 
groove feeding arrangement reduces the variation of viscosity 
due to temperature since the oil flow for a line source is 
three to four times that for other feeding arrangements at any 
given load.14 . The accuracy of the short bearing approximation 
is iliustrated in Appendix A where the error is negligible for 
L/D = 2, and is acceptable for L/D = $. 

The assumptions that undoubtedly need considerably more 
investigation are those that apply to the flow of the fluid 
through the cavity,under dynamic conditions, particulsrly 
where the striations join the downstream boundary. The steady 
state investigations(referenced in the Introduction) using 
transparent bearings have shown a finite and irregular distri- 
bution of striations, with the result that considerable mixing 
must occur at the downstream boundary, and thus in actuality 
the downstream boundary cannot be precisely defined mathema- 
tically. Also, because of the low pressure gradients immedi- 
ately surrounding the boundary, surface tension and wettability 
probably influence the shape of the boundary. 

The other assumptions that were made such as laminar flow, 
Newtonian fluid, etc., are generally accepted and will not be 
discussed in this paper. 
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CONCLUSIONS AND RECOMMENDATIONS 


This investigation has analyzed and programmed the un- 
steady motion of a circular, circumferentially grooved journal 
bearing. The results should prove useful in designing the 
particular type of bearing considered, while the analytical 
and computational methods employed can be modified for арр11- 
cation to other bearing configurations. The most important 
point that was brought out by this investigation was that ne- 
glecting the oil film between Ө = п and Ө = 2n will not 
yield good values for the spring and damping coefficients in 
lightly loaded high Speed bearings, where Ockvirk's theory 
that a given load establishes the attitude angle is invalid. 
Conversely, it would be expected that assuming a complete 
film in that region would also give poor results with the 
exception of hydrostatic bearings. 

Recommendations for future work in tne subject area of 
this investigation are listed below. 


l. Devise and construct an experimental apparatus to 
measure dynamic forces and load angles, with sufficient 
accuracy to establish the validity of the boundary 
conditions of this investigation as compared to 
Ockvirk's boundary conditions. 


2. Extend Richardson's theory and his program to obtain 
the spring and damping coefficients. If the differ- 
ences are reletively small, Richardson's analysis 
would be preferable due to the simplicity of his 
boundary conditions; in its general form it could 
more readily be modified for various particular anal- 


ysis such as half speed whirl, stability, vibrations, 
etc. 


3. For use with L/D>s, program Reynolds' three dimensional 
equation using the boundary conditions of contiuity. 
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APPENDICES 








APPENDIX А 


Accuracy of the Short Bearing Approximation 


Donaldson[ 7] obtained an exact solution to Reynolds! 
equation for a journal bearing using series solutions and ap- 
plying Sommerfeld's boundary conditions, but he assumed that 
the feed pressure was sufficiently high to suppress cavita- 
tion. Figure XVIII is a comparison of Donaldson's solution 
and the results of this investigation, using data for uncavi- 
tated bearings from Appendix C. Two length to diameter ratios 
are plotted as a function of eccentricity ratio versus the 
ratio of approximate to exact load capacity. 


Figure XVIII 


Approximate Load 
Exact Load 





Eccentricity Ratio 


Sl 





APPENDIX B 
Computational Procedure 


The analysis was programmed in Fortran II for use on the 
IBM 7094 digital computer. The program in its present form can 
be used either as a main program or as a subroutine by the in- 
sertion of the proper administrative and format statements. 

A sufficient number of comment statements is in tne program 

to permit someone who is familiar with the language to follow 
the logic, but some of the detailed computational methods need 
further explanation. 

It can be seen in Figure VIII that as the circumferential 
flow becomes mucn less than the axial flow tnat the downstream 
boundary conditions approach the upstream boundary conditions. 
and equation 22 applies throughout. If the flow is negative, 
mee, (UO - 20) is negative, then the positions of the upstream 
and downstream boundaries are reversed. То determine which 
condition prevailed, the following expression was used: 


circumferential flow А пе 


Velocity Ratío (VR) = а а |-Ano| «|&| 


This gave not oniy the direction of the ciroumferential fiow, 
but also its relative magnitude. If VR was positive, all 
iterations were done in the positive € direction, and if VR 
was negative, then all iterations were done in the negative 
Ө direction. 

The beginning of the upstream boundary was obtained by 
iterating in 3 degree intervals, starting at © = 7 and calcu- 
lating w by equation 22, until the following conditions were 


mes; 


ie wie L 
27-00 < w(1 -1) 


Ба 





If these conditions were not satisfied when the program had 
iterated through 2n, then a non-cavitated solution was obtained 
by means of definite integrals (see Appendix D). 

If cavitation was present, then 02 was located by iterat- 
ing in one degree intervals starting at 6, , again calculating 
м by equation 22, until w(I) > w(I-1) i.e., the sign of the 
Slope of the cavitation boundary changed. This established 6 . 
One degree intervals were used between 8, and 8, so that greater 
accuracy could be obtained in calculating h* later in the pro- 
gram. 

Then, starting again at Ө, ог ө; #5] whichever was even, 
the program was iterated in 6 degree intervals using equation 
22 to calculate the film width along the upstream boundary, and 
equation 27 to calculate the film width along the downstream 
boundary. Equation 27 was applied as follows: 


l. Assume that w(I) = w(I-1) 


2. Determine h* by iterating through the w's previously 
calculated between 6, and 6, , until w* > w(I). 
then interpolate between the stored values of 8% on 
either side of w(I) to obtain h*. 

3. Calculate the slope — by means of equation 27, using 
h*(I-1) and w(I-1). 


4. Then, w(I) => м(1-1) + = Ae 
However, because equation 27 becomes invalid, both mathemati- 
cally and analytically, as A.-0, and because the quantity 
(h-hiÀ*) becomes very small and practically impossible to calcu- 
late digitally with any accuracy as VR —~ 0, it was necessary 
to apply certain restrictions along the downstream boundary. 
These were: 


52 





l. If |VR 


| > i equation 27 applied. 
3217 IVA] =< 


425, equation 22 applied, 


The number .25 was obtained by trial and error; as VH app- 
roached this value, the downstream boundaries by equations 22 
and 27 agreed very closely (Figures XV and XVI have VR values 
Of .33 and -.33 respectively). At Lower values of VR , this 
program gave inconsistant results for the reasons mentioned. 

Por each © , the differential forces and coefficients 
were calculated by equations 28 and 28 and equations 37 through 
44, The integral of h*dz for use in calculating the friction- 
al forces in the cavitated zone was evaluated by Simpson's Rule 
femme 3 ordinates 


L 
E са 5 (h* - h*)(w* - E T W) 
Ww 

Then the differential forces and coefficients were integrat- 
ed from =m Бо ө = 3m, by means of the subroutine,which 
uses Simpson's Rule and 61 ordinates. For clarity, each compon- 
ent of each equation was calculated separately; thus 0С551 is 
Er zrirst term of IK... 

The remainder of the program is self explanatory. It 
should be noted though that FR is the radial force directed 
inward; this was changed for conventions sake. 
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N 
e 


ho 


r-? 
| 


FCRTRAN 
FORTRAN 


Priests ОР Лл CIRCUMPERENTIARLLY GROOVED JOURNAL BEARING 


EXULPSDISCTODNE BEARTNG LENGTH JOURNAL RADIUS RADIAL CLEARANCE, 
MUA  ЗКОГТАС SPpEEyRACIAL VELOCITY ECCENTRICITY- -RATIC, 
BIETE FEED PRESSURE, AWD VISCOSITY (EL,R,Cs0M,;PHIDOT, DE,EN, PO, U] 
ИМІТЗ ARE LD, IN ySEC RADIANS 


IPS RADIAL? TPANMGENTIALS,hORIZONTAL$,AND VERTICAL BEARING FORCES 
B GSC CCFERECCIESILSSAMD DAMPING COEFFICIENTS( FREY, ETOT,ANG CRE, 
E See Ps pb Oe tty CAXKG CYY s CY Xe GXKXV pt Xkre BYYV otk Ys ) 


ШЕШ ГОР Х2(99),М2199),Х(61),5(61), ОҒІ1(61),:О0Ғ121(61), ОҒІЗ(61); 
E LEIe[519,0F29*(6521);0-F22(61),;DF24(61); DCRR1(61) , DEORA2(61) ; DCSRSIUCÓOI Ms 
Bie 2161) ANSDZ (61), UGRKR3(6€1),0CRS4(61),00$84(61),DCSR3(61), 


Beet (61) yDF25>(61),W1 (124) 
BoMbine AND .KEDUCE INPUT 


GA=O0M—-2.0*PHi DOT 
Е-РО-15.0 
PI=3.1416 
ӘР — U.5*cL 

ШО - (0.>есі1)%%3 
BU-(UCTR)/U2.Cs8C x$ 3) 
Ср= (Оз ОМК 7) ИС 
SS=Du/C 

TT=CC=OM 


BS=0. 
20. 
pr=0. 
ШЕ-О. 
Cx=0. 
CS=0. 
CT=0. 
CU=0. 
СУ+0. 
CYzC. 
AHSOZ=C. 


E NITUDE ANC DIRECTION OF CIRCUMFERENTIAL FLOW 


ViX=GA/ (ABSF(GAJ)+ABSF(DE/(EN*C))) 
IF(VR) 2:0»,110,1:0 


meCcAiE ХІ FOR NEGATIVE FLOK BY ITERATING IN 3 DEG INTERVALS IN 
BeGAt{ VE THETA CIRECTICN,BEGINNING AT 5 PI 


ШІ 220 К=1, 
И (к=]25&)-.2 
A = K 
K=K 
X 
А 


124 
1‹,›300,999 


ISP 052395=(4.05235) : 
mi S=€=(1.0+oNeCOSF(K1)) 
>> 





BAA=COSF(XIL) 
CBB=SINF(X1) 


БОШАТТОМ 20 


PH= 175 <UsEL we 28 (-2 CEDE SAAA*GAe ENSCEBDO)/LCE (2.088 N8 AAA) ) 943 
Pelee OePii+P) 212,885,215 


ЕССАТТВМ 22 


212 wL(K)-SQRTF(-((CSUVL.O4*4ENSAAA) ) *3«P)/(3.C*U* ((72.0*DC* 
1AAA)J)+(GAFEN=CxBud)))) 


Mate) ISO) 220092 1 3 


eee (WLC A-G)-.c99* EL) 220,214,214 
Meet (W1i( K)—wWil(K-1)) 240,220,220 
8 ultK)20.5*EL 
pee CONTINUE 
BUS КО Fon NEGATIVE FLOW SY ITZRATINO IN CNE DEG  INTERVAVS UIN 
E AIVE THETA DIRECTION , GEGINNING AT Xl 
аи со 250 J=1,99 
J=J 
261 А=Ј 
X2(J)2X1-(A*.01745) 


СӘСАТІСМ 22 


BZ SORTE (UCERO. OFENFECSF1X21(3))))2=*3=P)/(3.0*U*((-2.0*DE E 
memes (X2(3)))+(GASEN*CESINF(X2(3)))3))) 
ШЕ(22-3) 250,250,245 
ELLOS (W2( 3) -w2(5-3)) 25C,251,251 
AA CUNT INUE 
Moe “e2=Vi26u-3) 
mee 2=X2( 5-3) 
а= 7-12 
шз6 1 -К/2 
Са ТО 153 


BECATE xi FOMEPOSÍ'IVE FLOW BY ITERATING IN 3 DEG INTERVALS IN 
Ш/ІТІУгЕ THEIA EIRECTICN,GEGINNING AT 3 PI 


ШОХ 0 120 K=1,124 
MmeetK-124) 181,300,999 

111 A=K 
K=K 
EE 3203634+(A2.05235) *2.0*PI 
ARH1IS=C*(1.0+EN=COSFIX1)) 
BAA-COSF(X1) 
EBB=SINF(X1) 


EQUATION 20 


Pri.18759*U£*EL &52*(-2.0*3DESAAA*CASENESC«EBBDI)/CCSUL.O*ENSAAA))* 932 
bet4.0*PmtP) 112,1153115 
56 





cQUATIGN 22 


р ~ 
к” 
rx 


i) КЕ = (Г [ОК О СТЕМЕДАА у УС ЕЗЕР) (3. 050+ 01-2 .0+05= 
1LAAA)*UGASENS C«BDE2)))) 


NDS, 205.113 
ШҮШГІГІНІ (К-:)-,499:Ғ51) 120,114,114 
ШЕ TE C WiC K)-WECK-1)) 140,120,120 
ШІ -Һ1(К)-0.5<сі 
OO CONTINUE 


Ee x2 FOR POSITIVE FLOW BY ITERATING IN ONE DEG INTERVALS IN 
ШІЗІТІУЕ ТНЕГА DIRECTION $, BEGINNING AT Xl 


140 150 J=1,99 


CO 

J= J 
141 А=)) 
Ха) <Х. + (А+.01 742) 
EQUATÍON 22 


ОО ЕБОКЕР СЕСЕ УСТЕМЕСОЅЕ(СХ2( 23) Уу 06 39р) (3. Ок уе ( (-2.0=DE ж 
ERSUSFIX2(J)JYTUGAE ENSC*S INF(X2(422 02020) 
г л—3) 150,650,145 
ESUIt(Ww2(J)-W2(J3-3)) 150,151,151 
ED Ceili INUE 
ЖІ WW2=W2 (5-3) 
XX2=X2(J-3) 
w= jJ-1 2 
I=K/2 
ШОО? К) 13.25 153,999 
152 I=I+1 


INCREMENT IN 6 DEGREE INTERVALS AROUND THE JOURNAL BEGINNING AT 
E AL EVER) Or AT Xx1+3(1F X- 0DD) 


ВО 170 N=1,62 

A=] 
X(1)=3.0#P1-.1047+(A%.1047) 
AAAZCOSF(X(1)) 
e£83=SINF(X(1)) 


ы?” 
Ал 
uz 


Meee LAT IVe FAGNITUDE CF CIRCUMFERENTIAL FLOW IS SMALL, 
ESUATION 22 APPLYS TRRUUGHCUT 


BECABSF(VR)-.25) 154,159,159 
Meee (VR) 2604,1604, 1604 


2604 IFUXX2-Xt 


I 4,154,166 
moo 1F(XX2-X(I 


JO 
AS ety lS 


EWUATION 20 


154 PH=.1875#UeEL e# 2% (-2.0#DE KA AA+GA® EN#C#BBB)/(C#(LIOTEN# AAA) 2 я3З 
PEAS ОжритР) 165, 167, 167 
ЭЕ 





ша SH y CHDUATIDMP 22) 


16» SUP)ESQATF(-((C(2.0*N&AAA) *832P)/(3.C08U& ((-22.09DE AAA) * (GA#EN 
1 #С+ В5В)) ) ) : 
BAHS=C* (2 .0+ENEAMA) 
XS=X(1) 
Eb TO 168 


CCHUNSTREAM W I$ COMPUTED FROM STATEMENT 166 THROUGH 1567 


Bere iF (VR) 11,15,15 
11 EN 51) 12,13,5999 
үр h(I)=w"(1+1) 
EC TO 20 
13 h(1)=Ww(1) 
ES TO 20 
ES IF(I-1) 999,16, 17 
JG К(СІ)-“1(61) 
ED TO 20 
17 h(1)=W([-1) 


ге РОК EQUATION 27 


20 EN 1564 4-1,3. 
BEC d-JJ) 36g 165,999 

36 ши )-и2(1)) 1664, 1661, 1661 

ШОТ ІГІ.)-1)99941662,1663 

BE XS-zx2(4J) 
E JU 1665 

EE -x2(J-210)*(0N(I)-u2(J—-1))«4(X2(J0)-X2(J0—-12)7(W2(2) -42(J- 12) 
CO TO 1665 

1664 CONTINUE 

moos AtiS=C#(1.0+ENeCOSF(XS)) 
AnzZCE(l.O*EN*AAA) 


EQUATION 27 
L&-(0.5*w(1)2(-2.0*DE*AAA*GA«EN*CSBBB)T((AHES*32P)/(6.08U&RV(10)0)) / 


(СА“(АН-АН5)) 
1667 W(I)-W(I)*TABSFUDW)4.1047 


b 


59 moen(l)-EL#0.5) 168,167,167 
Meee w(1)=0.D5*EL 
ШЕВ AA=(1.0+EN#A NA) #% 3 
BD=Ww(I)*=*3 
CO=(1.O+EN#AAA) #22 
RR=ÇGOQ* £ 2 
Уи+ (1. О+ЕМКААА) # 22 


BIFEFERENTIAL FUNCTIONS OF FILM WIDTH AND TRIGONOMETRIC FUNCTIONS 
ES CAVITATED ZONE 


ОЕ ЛЕ CARDS DOWN TO 1683 IF FRICTIONAL EFFECTS ARE TO GE 
ШОО ЕСЕР. 


1682 Ан502-0. 
BEN -A9SEL-AII)) 1683,1683, 1671 
ЖЕ! (ІПІ125-АН154( АН5-АН:5)/2.0 
ШИ-Т(АНТІ25/2)-1.0)/ Ек 





KX12-ACOSFDDU) 
Pee-sINF(ALZ) 


EWJATICN 20 


PBh=.182 75 *U+EL 


ше2е O:DE*SDCDD*GAZLCNECEEEE)/(CC8SCUT.O* ZN&DDD) ) 52 
PS Capra) 1672 2 


ШІ? ЕСЕ--ЕС 


fi 


EQUATION 22 


1673 W12-SQRTF(-((C* (1.0*ENSDOD)) s&235P )/(3.0*U*((72.085D0E* DDD) * ( GASEN 
meCe CCE) ) ) ) 


EQUATION 22 


WS=SORTF[=((L=( 1 .O+En=sCOSFI(XS)))=*32%P)/(3.0%L*((-2.0*DE+* 
LCOSF (XS) ) +(GAXENSCRSINF(XS))))) 


BSS-ABSF((AH.S-AHS)/2.0) 
AHSCZ <(555/3.0)# (1 (К) +4.Оки 12 + 15) 


1883 ССКАЗ( Г) +Ан5 ГЛ ВВБХААА /АА 
ПОКУ4 (Г) АНУСА + ААА / ЦС 
CCSSA( I) ZAHSUZ2«8BD/CG 
LCSR3(I)2AHSDZSAAA*22/AA 
CF15(I)2AHSDZ*BBO/QQ 
LF25(I)2AHSDZ*AAA/QQ 


ШІБЕЕКЕМГІЛІ ГОКСТІПНА5 OF FILM WIOIH AND TRIGONGME TIRIC FUNCTIONS 
meecOMPLETE FILM ZONE 


CFL1l(I)2(BB*sAAA«*2)/AA 

СЕ12 (Г) < (| ВВ“ ААА + ВСВ) / АА 
СЕ1З(Т) < 5 (Г) + ААА 

СЕТ 4 (1) <и (1) #СВЕ/ (1. С3ЕМ“АДА ) 
СЕ22 (1) - | ВБ+ ВВ « #2) / АА 
СР2З( Г) < (1) +БВВ 


СЕ? 4 (ГТ) <-и (1) ЕКААА/ (1, С+ Е Ме ААА ) 
ССКК1 (ГТ) - ВВ+ НОВ ДААД «+ г /КК 
ССКК (Г) <И( 1) + ДАА ЕВОВИСО 
ИООРІ( Т) = ВВЕ ГЕЕЕ 2 едда 
ССЅА2( ТІ) = М(І) АДА + +2 /0С 
Ше .49+Е-И(Г)) 1685,1685,42C 
420 м5174+6+ 1 
IF(M-360) 42; 422,421 
421 N=M-360 
422 Z72Z2-22.0*W(I)/cL 


БЕТТ ТТГ HERE IF THETA ANG FILM WIOTHCPER CENT) TN 
ШТТАТЕГ КЕСІСОМ IS LESIREC. 


ШО) ТЕ(МК, 168о, 1 658, 1 083 
1686 1-1-1 

ШЕСІ) 999 1687, : ГО 
1687 I=61 

EU TO 170 
S38 I=1+1 

IF(I-62) 170,1689,999 
1689 I=} 


p 





Eu qe ЛО ЖОГ ЕЛЕУ ЛЛ ТА АХ DIFFERENTIAL TRIGONOMETRIC 
СЕТЕ ТОИ 05 ИС SIMPSONS RULE 


ЕД =SIMP(DFI,) 


EC = SIMP(DF12) 
mie = SIVPIUELEE) 
PG = SIMP(CF.2) 
ËH = SIMP (Uhm) 


ЕК-СІМР(СОСКЕ ) 
E -SIMP(DCSRS 9) 


ЖЕНЕ ТЕАТ ТӘТЕ FRICTIONAL EFFECTS ARE- TO BE NEGLECTED 


ECSIMPLDCSRZ) 
ESESIMP(DCRRAJ 
EE = SIMP(DF14) 
Bro = SIMPLDF24) 
R=SIMPCDCRR3) 
CS=SIMP(DCRSS) 

Br=SIMP(CCSS+) 
EUSIMPIDCSR2J 
ШИВОТМР(ОЕ:5) 

ШСІМРІСЕ25) 


COMPONENTS OF EQUATIONS 28,29;AND 37-44 


FLll=-2.0#CC#DEeBA 
F12=CC*GA*EN *C# BC 
F13=0.5%P#ReUD 
Fl4--DD*8E 

F152-SS*DV 
F21-2-2.0*CC*DE« BC 

Е22 СС #БА#ЕМ« Си ВС 
F23=0.5:P+*R*BH 

Е24=00 +80 

F25=SS=DY 
CRR1=TT*(8C-3.0+EN*BR) 
CRR2=SS#BS 
CRR3=2.0*SS*1R/C 
BRR=-2.0#CCEBA 
CRS1=-TT*BG 
CRS2=-P+R*BH/2.0x Cr ER 
CRS3=-SS*BP/EN 
CRS4=-SS*DS/CHEN 
BRS=-2.0*CC*3C 

oS 1=TT#0C 
С552-Р«ДаВб/2.0%С«ЕМ 
CSS32-SS :BE/.N 
CSS4=-SS*DI/C#EN 
RBSS--2.0*CC*8G 
CSR1-TT«4(80-5.0«EN«BV) 
CSR22-SS «BH 
CS$R32-2.0«S$$«DU/C 
BSR--2.0«CC«8C 


ED TO. 350 





Bn cc On stone CAMPING CONSTANIS IF NO CAVPTATION 


BOO —ERINT 9i? 
Alz-PI*EN/(1.0-CN&**2)*52.5 
AOA age lean со СЕМЕ 2660.5 ) ) 
АС = teof/().0-CNee2 ув | 5 
FilzO. 
F12-2-0.0 
БІЗ--0.0 
Ғ14-0. 
ЕІ5-0. 
F21=0. 
F22-CC*GA&*&GG«*EN«2.07AC&C 
F23=0.0 
F24=0. 
F25=0. 
CRR1=0. 
CRR2=0. 
CRR320.0 
ERR=0.0 
Са51=-бС»Т +2” .0*АС 
С452-0.0 
и: 53-- (554 ЕЕ/ЕМ ) #2 ОЗАО 
CRS4=0.0 
ER S-0.0 
55 1=0.0 
ЕЗ52-0.0 
С553-0.0 
CSS4=0.0 
pSS--GG*4.0£*CC* AC 
ESRI-GO*TT&( AC S3.OSENSAI) 
ErR2=0.0 
ESR3=0.0 
ESR=0.0 


NNF R=-=2.0+(F11+F12+F13+F14+F15)/([(7.5r+EL=R) 
EE Он ЕО te eet, 35rF24t4$r25)/(7/.5vVELSR) 
FETOT=SORTF(IFi:=2+FR=x=2) 

FHIZATANF(FT/FR) 

ANG=57.3*PHI 

en = 27.08 ERRIEERFZZFLRR3 JeC/( 7.5*EL eR) 

ENS --2.05(CRSIDECREZECRSS3TCRSA4 )*C/(T7T.5«*EL*R ) 
ша - 2.0 (СЗ 16 Е СЪ5З 554 )eC/( T. 5*EL «R ) 
Mere —-2.0*(CSRIFCSREFCSKS JRC / (7.52 EL ee ) 
PRR=-2 .O#BRREOM#C/( 7.58 EL #R) 

Boo=—2 «Or OMeC/( 7. S#ELeR ) 
ERS=2.0*BRS*OM*C/[7.9*EL*R) 

BAA-COSF(PHI)s5Z 

EDB-SINFC(PHIJ)s52 

ECC-COSF(PHI)*SINFUPRBI) 

ШАХ -СККедАд+ ( СЗК-СЕ5) + СССТ+С55 + БВБ 

ES CSNEBSBt(CRS-CSROIFCUCCTCSSt*AAA 

E YSCRSSAAAT(CRR-CSS)SCCCTCSR&eBBB 

NO URSSDODÓSJCULCSSCCRROSCCOCSECSRS*AAA 
EXX=BRAREAAA+?.0*BRS=:CCC+BSS+BBE 

EN Y=BRRF+BBB-=2.0*BRS=CCCFBSSRARA 
EXY=BRR=*CCC-=BRSH+AAA-ESS*CCC+BRS:BBB 


993) CALL EXIT 
END 
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FORTRAN 
SIMPSONS RULBRFGR 6? CREBLNATES,-MITH 6 OEGREE INTERVALS 


FUNCTION ЗІРЕЖҮ? 
Ша Е 45 LON Y( 1) 
CuD=0. 

EVEN=0. 

J СОО 1=2,60,2 
EVEN-EVEN-*Y(I)J 
СО 810 1=3,59,2 
СООЕОБС+У (Г) 
SIMP=(.1047T1/3.0)#(Y(1)+%#.O#EVEN+2.O#ODD+Y(61)) 
RETURN 

END 

CATA 


62 








APPENDIX C 
Tabulated Values of Bearing Load, Attitude Angle, 
And Spring and Damping Constants 


Tables 4 through 45 contain values of dimensionless load, 
load angle, and dimensionless spring and damping constants as 
tion of eccentricity ratio, pressure ratio, and dimension- 
Mees squeeze. The eccentricity is in increments of 0.1 from 
ЕКЕ to 0.7, the latter being chosen as the upper linit in 
Order to remain within the limits of the short bearing approx- 
imation, the assumption of constant viscosity, and the assump- 
mos of no journal or bearing distortion. The pressure ratios 
are 1:5, 2.0, 3.0, 4.0, 5.0, and 6.0, which adequately cover 
the range of most bearing applications. The choice for varia- 
tion of squeeze was made as follows: from Shaw and Macks 


R 1 
small bearings: = £ —— 
soot .003 
R i 
large bearings: — > ——— 
š a, 77770005 


Now let е = e Vary from 2 to l. Then R varies from 2L to 4L 
and: 





г 
file 1 |" _ 5 
“мім 
E s [ E 
E 277 7 
ES. | 0005 | 


Now assume a viscosity range from 10° to10° lb-sec/in? and 
a speed range from 100 to 2000 rad/sec. Then: 


dE 


АЕ + = 2)2 
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# = 077109912 8)(10* ) . 185 


МАМ 


А T (10 ) (2000) (1.6) (10 ) В 51590 
MAX 15 


However, no bearing applications at the upper limit could be 
conceived and so it was decided to vary A* from „1 to 8000 
in 30 steps that would permit easy interpolation. It should 
be noted that generally the first few entries in each table 
are uncavitateda bearings, wnich can be recognized by virtue of 
the fact that f = 90 

The tabulations are in terms of the entire bearing length 
(2L) and are converted to dimensional values as follows: 





F= F p LR 

Y шш NT С. 

E nm ANE DIR 
Ue 

Ва m= BNM р, LE 


For end fed bearings of total length L, divide the above 
values by 2. 
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APPENDIX D 
Table of Integrals 


The following definite integrals are used in evaluating 
equations 28 and 29 and equations 37 through 44 for the non- 
cavitated case. 
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